Intermolecular Adhesion in Conjugated Polymers 
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Conjugated polymers are observed to aggregate in solution. To account for this observation we 
propose a inter-chain binding mechanism based on the intermolecular tunneling of the delocalized ir- 
electrons occurring at points where the polymers cross. This tunneling mechanism predicts specific 
bound structures of chain that depend on whether they are semiconducting or metallic. Semicon- 
ducting chains should form polyacene-like states exhibiting binding at every other site, while (doped) 
metallic chains can bind at each site. We also show that solitons co-localize with the intermolecular 
binding sites thereby strengthening the binding effect and investigate the conformational statistics 
of the resulting bimolecular aggregates. 
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PACS numbers: 82.35.Cd,83.80.Rd,61.25.Hg 

Conjugated polymers are of great fundamental in- 
terest in that they are exotic, low-dimensional (semi-) 
conductors^ Q- More recently the use of such poly- 
mers has been explored in a variety of applications 
ranging from optical @, 3 and electronic Q devices to 
biosensors (3- While exploring their use as highly sen- 
sitive biosensors, researchers |7j discovered their surpris- 
ing tendency to aggregate into multi-molecular filaments. 
This aggregation is even more surprising in that it occurs 
in a system of strongly charged conjugated polymers dis- 
solved in low ionic strength aqueous solutions. 

Motivated by this observation, we investigate more 
generally a potential aggregation mechanism that should 
play a role in any 7r-conjugated system. In this letter 
we propose that the appearance of intermolecular elec- 
tronic states made possible by 7r-electron wavefunction 
overlap where two conjugated polymers cross leads to a 
large enough reduction of the electronic energy of the two 
polymer system to cause intermolecular binding even in 
dilute solution. The electronic intermolecular interaction 
draws localized states out of the continuum of extended 
states on each chain in a manner analogous to the for- 
mation of traditional bonding/anti-bonding orbitals. We 
examine the effect of this bonding mechanism in both 
semiconducting (SC) and doped metallic (M) polymers. 
In the case of the SC polymers we also examine the in- 
teraction of solitonic excitations on the chains with the 
binding sites. Finally we use the results of these analyses 
to predict the aggregated fraction of chains in solution 
and the change in the polymer conformational statistics 
that may be measured via scattering. 

We find that the inter-chain electronic tunneling leads 
to a significant change in the electronic states far from 
the Fermi level. In SC polymers, this change makes en- 
ergetically favorable a polyaclene-like conformation (see 
Fig. [21) with binding at every other site. In addition, there 
is a strong interaction between these inter-chain binding 
sites and solitonic excitations on the chain. Uncharged 



solitons are strongly attracted to the interchain binding 
points and, at room temperature, co-localize there fur- 
ther strengthening the inter-chain attractive interaction. 
Charged solitons, on the other hand, interact weakly 
with the tunneling sites at inter-chain crossing points. 
In M chains, however, the inter-chain electronic interac- 
tion favors the formation of a parallel polymer configu- 
ration where there is electronic tunneling at every site 
(see Fig. [2J and there are no solitonic excitation on the 
chains. 

In either case the presence of intermolecular bonds 
changes the statistical properties of the conformations 
of the bound pairs primarily by modifying their effective 
persistence length. As the binding strength increases, the 
dominant deformation mode of the polymers goes from 
bond rotations in the free chains to localized bending at 
defects in the intermolecular bonding pattern for weakly 
bound chains, and finally to the elastic bending of stiff 
polymer bundles. 

To consider the most simple description of conjugated 
polymers, we study a tight-binding model of noninteract- 
ing electrons based on the SSH Hamiltonianp] for each 
polymer chain. At a specific set of sites a (i.e. the cross- 
ing points between chains) we allow inter-chain tunneling 
with the overlap integral t' ~ 0.2eV ||. Thus, we write 
for the single particle Hamiltonian 

H = ^[-tt*+i{\t,n){t+l,n\ + \t+l,n)(e,n\)] 

e.n 

- i'£[M)(a,2| + |a,2)(a,l|]. (1) 

{a} 

Here the \t, n) represents an electron on the i^ 1 tight- 
binding site of the rfi 1 chain where the sum £ = 1, . . . N 
extends over all sites and ./V is the polymerization in- 
dex and n = 1,2 labels the two chains. In the case of 
the undoped, homogeneously dimerized chain t^+i = 
h + We expect t ~ 2.5eV and A ~ O.5eV0. 
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For homogeneous the M chains (i.e. without the Peierls 
distortion) A = 0. 




-3-2-10 1 2 3 

Re(2ka) 



FIG. 1: (color online) Fhe spectrum of electron states result- 
ing from electronic intermolecular tunneling at the crossing 
point. The bound states are represented by the (red) filled 
circles at Re(fc) = 0, ir. The schematic figure in the center 
represents such a single tunneling site between two crossing 
polymers. 

We first consider a single crossing point between two 
homogeneous, infinite, undoped chains using a transfer 
matrix approach 9] . The Hamiltonian of the system com- 
mutes with the chain interchange operator so that the 
energy eigenstates have definite parity with respect to 
chain interchange. We find four bound states, two in 
each of the symmetric and antisymmetric sectors of the 
state space, at energies 



(2) 

where we have defined: ti^ = fo ± |. These bound 
states are represented as red circles on the band struc- 
ture diag ram shown in Fig. ^ There are two ultra- 
band [lOj states above and below the valance and conduc- 
tion bands. These states are split from the band edges 
by 0(t' 2 /to) and their wavefunction decays over O(to/t') 
sites from the tunneling site. Due to this localization of 
the bound state wavefunction, two tunneling sites sep- 
arated by more than this distance interact weakly. In 
addition, there are two states that appear symmetrically 
shifted by C(t /2 A/tg) into the gap. They have a real 
part of their wave number at the edge of the zone and 
localization length Oiat^jt' A). The effect of the polymer 
ends on all of these localized states is insignificant unless 
the tunneling site is within a localization distance of the 
chain ends. For high molecular weight chains, we may 
ignore end effects. 

The dominant part of the reduction of the system's 
electronic energy is due to the occupation of the lower 
ultraband state; the contribution of the gap state can 
be neglected. Similarly, the shifts in the extended states 
of the chains are small for high molecular weight chains 
since these shifts vanish as 0(N^ 1 ). Previous calcula- 
tions on the M chains 11] demonstrate that the ultra- 
band state is displaced by 2t[ x /l + t' 2 /(4t 2 ) - 1] below 
the bottom edge of the band, while, once again, the effect 
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FIG. 2: (color online) The energy of the bound undoped 
chains as a function of the fraction p of (equally spaced) bind- 
ing sites for t' = 0.25eV. The insets show (left to right) two 
free chains (p — 0) , and the expected forms of the maximally 
bound SC chains at p = 1/2 and M (doped) chains at p — 1. 
The former structure is similar to polyacene except the bind- 
ing results from the hybridization of the n z orbitals. 

of the shift in the continuum states and the gap state are 
irrelevant. In the case of both the M and SC chains, we 
note that the appearance of a single localized state cen- 
tered at the crossing point of two polymers reduces the 
electronic energy of the chains by ~ lOOmeV. 

We now consider the most energetically favorable con- 
figuration of the two chains in both the M and SC cases. 
Here we demonstrate a distinct difference between these 
two cases. Previous work on the M polymers has shown 
that the strongest interchain binding in the absence of 
chain electrostatic repulsion occurs for the parallel chain 
configuration as shown on the right side of Fig. El In this 
case the individual localized ultraband states below the 
band edge hybridize into a new filled band. As this result 
seems intuitively reasonable, it is surprising to find that 
in the case of the SC chains the parallel configuration of 
the two chains results in no binding energy at all. 

One can understand this absence of intermolecular 
binding in the parallel configuration as follows. In this 
configuration the inter-chain interaction symmetrically 
splits the filled valance band of the SC polymer. Since 
the band remains well-separated from the similarly split 
conduction band there is no net energy reduction upon 
splitting. For the case of strong enough inter-chain tun- 
neling (t' ~ A) so that the valance and conduction bands 
do not remain well separated, the Peierls distortion of the 
chains breaks down and returns the system to the case of 
M polymers 8]. Instead we find that the energy minima 
for SC chains is the polyacene- like structure shown in the 
middle cartoon in Fig. [21 

To explore this further we compute the electronic con- 
tribution to the binding energy by direct diagonalization 
of the two-chain Hamiltonian with 200 tight binding sites 
on each; the results are shown in Fig. [21 This figure is 
symmetric about the p — 1/2 tunneling site density that 
is the minimum energy configuration. The symmetry be- 
tween binding site densities p and 1 — p can be under- 
stood as follows. Removing a single binding site from 
the fully linked chains can be represented by adding an 
inter-chain tunneling term to Eq. ^with t' — > — t'. Thus 
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FIG. 3: The upper figure shows the interaction energy of 
two solitons one on each chain separated by s' lattice sites. 
One soliton is centered at the binding site. The lower figure 
shows the interaction energy of a pair of symmetrically placed 
solitons displaced from the binding site by s. 



a single non-tunneling defect in the fully linked chains 
generates the same set of localized states as does an iso- 
lated tunneling site on the unlinked chain. From Eq. [3 
both localized states lower the electronic free energy of 
the system by the same amount. We extend this argu- 
ment to multiple non-tunneling defects in the otherwise 
fully linked chains by noting that the associated local- 
ized states simply hybridize to a degree that depends on 
their spatial separation. At any finite density of regu- 
larly spaced tunneling (non-tunneling) defects the local- 
ized states form a narrow "impurity" band, but the sym- 
metry between the two cases is unchanged. The effect 
of quenched spatially randomness in the distribution of 
the tunneling (non-tunneling) sites on the energy of the 
system is minimal 

In SC chains the inter-chain binding interaction is com- 
plicated by disorder in the dimerization pattern that 
takes the form of solitons. These defects may arise due 
to either end effects, cis- insertions in the trans- chain, 
or light electron/hole dopingflUT^. We now explore the 
interaction between these solitons and the binding sites 
on the two chains by considering two chains linked at one 
tunneling site and examining the energy spectrum of the 
system as function of the relative positions of the solitons 
from the binding site. 

First, we posit that each chain contains a single soli- 
ton and vary the position of a the inter-chain tunneling 
site with respect to the positions of the soliton. See the 
inset of Fig. [^Jd for schematic diagram. In moving the 
tunneling site, we require the tunneling site to remain on 
equivalent sites on the two chains to preserve the simulta- 
neous diagonalization of the Hamiltonian and the chain 




FIG. 4: (color online) The mean length (L z i p ) of the bound 
regions versus the binding E\> energy between tunneling sites 
on the two chains. The inset shows the fraction ji of bound 
chains as a function of Et in dilute solution. 

interconversion operator. Initially, we approximate the 
soliton as a single-site shift in the dimerization. We con- 
sider a more physical model of a spatially extended soli- 
ton below, but this simplified model captures the basic 
effect of the solitons in an analytically tractable manner. 

The single site soliton introduces new localized states 
at energies E given by the roots of the polynomial, 

= E 3 - E{A[t 2 + A 2 ]+t' 2 )±8t'A, (3) 

where the upper (lower) sign applies to the chain anti- 
symmetrized (symmetrized) states. For large E we find 
bound states at E ~ (4[tg+A 2 ]+i' 2 ) 1 / 2 similar to the pre- 
viously discussed ultraband states. For small E, however, 
we find a pair of gap states with energy E ~ ±i' A /to that 
are symmetrically positioned about the mid-gap of the 
SC chain. These are remnants of the usual soliton asso- 
ciated mid-gap states of the two chains split by the tun- 
neling matrix element t' . The mid-gap state associated 
with an uncharged soliton on the chain is singly occupied 
so that the splitting of these states allows for the double 
occupancy of the lower state resulting in a reduction of 
the system's energy upon the creation of the inter-chain 
tunneling site. The localization of an uncharged soliton 
on each chain at their crossing point further lowers their 
electronic energy by A/t' ~ 2.5 enhancing the binding 
energy at that site. In contrast, for the case of charged 
solitons (as would occur on lightly doped, SC polymers) 
both mid-gap states on each chain are filled there is no 
further energy reduction due to the co-localization of the 
charged solitons with the binding site. 

To extend this result to more physical solitons we 
put one SSH soliton of width la at the center of each 
chain and numerically calculate the binding energy for 
the TV = 199 chains. The interaction energy between the 
SSH solitons and the tunneling site is shown in Fig. [3] 
We find that the solitons are attracted to the inter-chain 
tunneling site. The width of this potential well for soli- 
tons is broadened once one soliton reaches the binding 
site. The interaction potential between the solitons and 
the binding site is oscillatory because the soliton associ- 
ated mid-gap states have nodes on alternating sites. 

The formation of these bimolecular assemblies will ef- 
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feet the photophysics of the system[13| as well as the 
conformational statistics of the polymers in solution. As 
the site binding energy Eb is increased there is an abrupt 
condensation of the polymers in dilute solution into bi- 
molecular filaments. See the inset of Fig. 0] We do 
not consider aggregation into larger supramolecular as- 
semblies. The persistence length of the bimolecular fil- 
aments has three distinct regimes as a function of Eb- 
In the limit of a vanishing Eb, conformations of the free 
chains may be approximated as a freely jointed chain 
with a Kuhn length of l m i n ~ 10a^^ where a is the 
length of a unit cell. For Eb ^S> k^T the bimolecular 
filaments saturate their available binding sites, and the 
intermolecular bonds hinder bond rotations so that the 
softest deformation mode involves bending in the conju- 
gation plane. Using known bond bending and stretching 
constants [H Hq| . we estimate that the persistence length 
is Lax ~ 100a. Between these extremes, where defects in 
the intermolecular binding are typically separated by less 
than £ max , the chain tangents decorrelate primarily due 
to the essentially free rotations at those defects. Thus 
the bimolecular filaments will have a continuously vary- 
ing persistence length £p, £ m - m < £p < £ max and their 
effective persistence length is determined by the mean 
length of the fully bound chain segments (L z i p ). 

We solve for this mean length (£ z ip) using a modified 
Poland-Scheraga (PS) model 16]. The partition function 
for the partially bound pair of chains each containing N 
sites may be written as 

p {ii 1) A 2 \h} 7r=0 

where u(i^\ i^) is the Boltzmann weight of the 7r th un- 
bound region consisting of i„ + tight binding sites 
where the superscript indexes the two chains. v(j n ) is 
the corresponding Boltzmann weight associated with a 
bound region of j w consecutive sites. The sum in Eq. 0] 
is constrained by the fixed total number of tight binding 



sites. We relax this constraint using the Grand Canoni- 
cal ensemble and, from the PS technique we compute the 
free energy F — — ksTN \n(x{) where x\ is the root of 

xv-v 2 {l-a) ^ \x) \2J ■ 1 ' 

Here v = exp(— Eb/k-oT) is the Boltzmann weight of a 
single binding site, while u ~ exp(a/£ p ) « 1 is the expo- 
nentiated entropy of a segment of unbound chain. The 
exponent c ~ 2.1 reflects the loss of entropy in these un- 
bound regions due to self avoidance and the loop closure 
condition [13 . Finally, -k B T\n a is the domain wall en- 
ergy between bound and unbound regions accounting for 
local chain bending energy. The fraction of binding sites 
is given by 9 — dhxxx/dhxv and the mean length of 
bound regions (L z i p ) may be computed from the ratio of 
8, to the fraction of sites at the boundary of bound and 
unbound regions, 6b- The results are shown in Fig. 0] 
For E b < 0.4fc B T, (L zip ) < £ max . Thus, given the best 
estimate of t' and for N > 10 3 , chains will aggregate in di- 
lute solution (10 _5 M) into this intermediate state where 
the persistence length is continuously controlled by Eb- 
Measurements of the persistence length may test the es- 
timates of t' arising from quantum chemical calculations, 
however, there are other contributions to Eb including the 
electrostatic repulsion between e.g. the highly charged, 
water soluble polymers. More direct tests of this the- 
ory necessitate measurements on more chemically simple 
^-conjugated systems such as polyacetylene. 

Understanding the persistence length of bimolecular 
filaments may play a role in controlling the penetration 
of such polymers into microporous or zeolitic structures. 
Future work will explore conjugated gels formed at higher 
polymer concentrations. 
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